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Abstract 

P-Glycoprotein, a member of the ATP-binding cassette (ABC) superfamily, is a multidrug transporter responsible for cellular 
efflux of hundreds of structurally unrelated compounds, including natural products, many clinically used drugs and anti- 
cancer agents. Expression of P-glycoprotein has been linked to multidrug resistance in human cancers. ABC transporters are 
driven by ATP hydrolysis at their two cytoplasmic nucleotide-binding domains, which interact to form a closed ATP-bound 
sandwich dimer. Intimate knowledge of the catalytic cycle of these proteins is clearly essential for understanding their 
mechanism of action. P-Glycoprotein has been proposed to hydrolyse ATP by an alternating mechanism, for which there is 
substantial experimental evidence, including inhibition of catalytic activity by trapping of ortho-vanadate at one nucleotide- 
binding domain, and the observation of an asymmetric occluded state. Despite many studies of P-glycoprotein ATPase 
activity over the past 20 years, no comprehensive kinetic analysis has yet been carried out, and some puzzling features of its 
behaviour remain unexplained. In this work, we have built several progressively more complex kinetic models, and then 
carried out simulations and detailed analysis, to test the validity of the proposed reaction pathway employed by P- 
glycoprotein for ATP hydrolysis. To establish kinetic parameters for the catalytic cycle, we made use of the large amount of 
published data on ATP hydrolysis by hamster P-glycoprotein, both purified and in membrane vesicles. The proposed kinetic 
scheme{s) include a high affinity priming reaction for binding of the first ATP molecule, and an independent pathway for 
ADP binding outside the main catalytic cycle. They can reproduce to varying degrees the observed behavior of the protein's 
ATPase activity and its inhibition by orfho-vanadate. The results provide new insights into the mode of action of P- 
glycoprotein, and some hypotheses about the nature of the occluded nucleotide-bound state. 
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Introduction 

The multidrug transporter P-giycoprotein (Pgp, ABCBl) is a 
plasma membrane protein belonging to the ABC superfamily 
which couples the efflux of a wide variety of chemically and 
structurally difiFerent compounds to the hydrolysis of ATP [1]. 
Commonly used chemotherapy drugs are transported by Pgp, and 
its overexpression in tumour cells is linked to the multidrug 
resistant (MDR) phenotype that many human cancers present in 
the clinic [2,3]. Following the first report of Pgp ATPase activity 
[4], studies characterizing ATP hydrolysis were conducted in the 
early 1 990s with plasma membrane preparations from MDR cell- 
lines [5], partially purified [6,7] or purified detergent-solubUized 
Pgp [8,9], and reconstituted Pgp [10,11]. Since then, the catalytic 
cycle of the enzyme, its coupling to drug transport, and its 
inhibition by oriAo-vanadate (Vi) have been studied by several 
research groups [12,13]. 

In 1995, Senior's group published a minimal reaction pathway 
for hydrolysis of one molecule of ATP by Pgp, and Vj-induced 
inhibition of its catalytic activity [14]. The protein possesses two 
consensus sequences for ATP binding, located within the two 



nucleotide binding domains (NBDl and NBD2) in the highly 
homologous halves. In support of the proposed scheme for the 
catalytic reaction, it was demonstrated that both NBD 1 and NBD2 
are capable of binding and hydrolysing ATP [14-16]. Thus, the 
minimal reaction scheme presented for the hydrolysis of ATP and 
trapping by Vi corresponds to the catalytic activity carried out 
independentiy by each half-molecule. Consequently, the apparent 
single K,„ observed for ATP hydrolysis [5,9], and the apparent 
single Kj reported for binding of nucleotides and nucleotide 
analogs observed by fluorescence and EPR spectroscopy [17-20], 
suggest that NBDl and NBD2 are similar in their binding and 
kinetic properties in regard to the hydrolysis of ATP. 

It is now generally accepted that the two NBDs of ABC proteins 
must interact to form a sandwich dimer for the normal functioning 
of these proteins, and such cooperativity has been shown for Pgp 
[21]. Thus, inactivation of one of the catalytic sites by either 
mutation [22] or chemical modification [15], or the formation of a 
non-covalent long-lived complex with V; trapped at a single NBD 
[14,23], is enough to completely abolish the ATPase activity of the 
enzyme. As result, steady-state catalysis takes place only when both 
half-molecules are intact. In addition, stimulation of the basal ATP 
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activity by drug substrates is displayed only when the full-length 
transporter is expressed, or both half-molecules are co-expressed 
[24]. 

Based on a minimal reaction scheme, Senior and co-workers 
[25] were the first to postulate a model for coupling between the 
catalytic activities of the two NBDs, under the name Alternating 
Catalytic Mechanism. In this model, the hydrolytic reactions of each 
half-cycle, and the accompanying protein conformational changes, 
alternate to carry out the transport of a drug molecule. The 
catalytic activity at NBDl containing a bound ATP molecule is 
triggered when a second ATP molecule binds to NBD2, and vice 
versa. Thus, the reaction progresses in an alternating sequence of 
ATP binding and hydrolysis in the complementary half-molecules. 
In this model, the authors hypothesized that the transport of the 
drug is coupled to the relaxation of the protein from a high 
chemical potential state that is generated by the hydrolytic step. In 
one turnover of this cycle, two molecules of drug are transported 
and two molecules of ATP are consumed. This basic mechanism 
for the catalytic and transport cycle of the Pgp is currently widely 
accepted, with the addition of further adaptations based on 
structural and energetic considerations. 

The Alternating Catalytic Mechanism suggests that asymmetry in the 
two halves of Pgp must be maintained throughout the catalytic 
cycle, in order to retain the memory of which NBD recently 
hydrolysed ATP (for a detailed discussion of the role of asymmetry 
in ABC protein function, see [26]). Using Pgp carrying a mutation 
in an essential catalytic residue in both NBDs (E552A/E1197A), 
Tombline et al. demonstrated the existence of a stable asymmetric 
nucleotide-bound Pgp species [27]. After gel filtration chromatog- 
raphy the protein retained one molecule of ATP, which was bound 
with 50-fold higher affmity {Kj— 9 |J,M) compared to wild-type Pgp 
(A^=0.5mM). This asymmetric species with a single tightly- 
bound ATP molecule was referred to as the occluded state. Sauna 
and co-workers reported the occlusion by wild-type Pgp of a single 
molecule of ATPyS, a very slowly hydrolysable nucleotide analog 
[28]. More recently, our laboratory reported that wild-type Pgp 
binds ATPyS to form an asymmetric species with one tightly 
bound (occluded) nucleotide {Kj=& |a,M) and one loosely-bound 
nucleotide {Kj=QJ mM) [29]. The asymmetric intermediate was 
proposed to exist transiently during the catalytic cycle, with the 
occluded nucleotide normally undergoing rapid hydrolysis [27,29] . 

In spite of the great advances in understanding the ABC 
superfamUy that have taken place over the last 15 years, no 
comprehensive kinetic analysis has been carried out to date. The 
compact mode of catalysis proposed by Senior's group in 1995 has 
been used to establish possible transport mechanisms (e.g. it was 
used in the Sequential Mechanism proposed by Sauna and Ambudkar 
[30]), regardless of the fact that several puzzling experimental 
observations have only been described superficially, and no 
satisfactory explanation has yet been proposed for them. These 
previously ignored observations could possibly be key pieces of 
information in the development of a comprehensive kinetic model 
for the catalytic cycle of Pgp. In this work, we built several 
progressively more complex kinetic models, and then carried out 
simulations and detailed analysis to test their vahdity in the 
proposed reaction pathway for the Pgp-mediated hydrolysis of 
ATP and its inhibition by V;. To establish kinetic parameters for 
the catalytic cycle, we made use of the large amount of 
accumulated data on verapamU-stimulated ATP hydrolysis by 
hamster Pgp, both purified and in membrane vesicles. We show 
that the proposed kinetic scheme(s), which include additional steps, 
can reproduce to varying degrees the observed behavior of the 
protein's ATPase activity and its inhibition by Vi. The results 



provide new insights into the mode of action of Pgp, and some 
hypotheses about the nature of the occluded state. 

Methods 

Construction of the Elemental Cycle Kinetic Model 

The basic kinetic cycle consists of an adaptation of the one 
originally proposed by Senior and co-workers in 1995 [14,25], 
here called the Elemental (Catalytic) Cycle (Figure 1). In it, a single 
reaction for binding and hydrolysis of MgATP is followed by 
sequential release of the products Pi and MgADP. For brevity, 
MgATP and MgADP will henceforth be referred to simply as ATP 
and ADP. Inhibition by Vi is achieved by formation of a long-lived 
complex, with ADP-V; trapped in one catalytic site, by a single 
step. This complex is thought to resemble structurally the normal 
transition state conformation formed with P;. In our implemen- 
tation, all the reaction steps were considered reversible except for 
the ATP hydrolytic step, which is irreversible [31]. 

The strategy used to test the validity of the reaction scheme in 
the Figure 1 consisted of the development of a kinetic model 
based on the rate law for each reaction. For modeling the scheme 
in the Elemental Cycle, the system was considered as a reaction 
medium without compartmentalization (which is the case for in 
vitro assays using solubHized enzyme or membrane vesicles) 
composed of the enzyme (E, Pgp) and one or more ligands 
(substrate, ATP; inhibitor, Vi; and products, ADP and Pj), for a 
total of n = 5 enzymatic intermediates. From the scheme in 
Figure 1, we formulated a set of ordinary differential equations 
for the rate of change of the concentration of n- 1 intermediates; 
considering each reaction as an elementary mechanistic step. The 
reaction system was defmed by a vector C„ of initial reactant 
concentrations, the total concentration of enzyme, [£];, and a 
vector k constituted of unimolecular and bimolecular rate 
constants, according to each unidirectional reaction, given by 

c,=<XATPUADPUPunoy 

with <[^rP]„,[^Z)P]„,[P,]„and[K,]„. The two following 

biochemical variables were solved either symbolically or numer- 
ically: 




Figure 1. The Elemental Catalytic Cycle of Pgp and Vi-induced 

inhibition. This scheme for the basic catalytic reaction for ATP 
hydrolysis by Pgp is adapted from Urbatsch et al. [14] E: Pgp. 
doi:10.1371/journal.pone.0098804.g001 
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(i) Turnover rate (in s ), v 



v = v(k,Co,[i?],) = 



. k2[E^ 



[Eh 



k.Co 



(ii) Fraction of trapped enzyme (adimensional), T 



T = T(k,Co,[E],)-- 



[E\, 



k,Co 



For modeling the other reaction schemes derived from the 
Elemental Cycle, in what are called extensions of the Alternating Cycle, 
additional differential equations were included to account for the 
new intermediates. In this regard, for the Partial-Extended Alternating 
Cycle (Figure 2, including blue reactions) the variables v and Tare 
defined by the expressions 



= v(k,Co,[P],) = 



r=T(k,Co,[Pl,) = 



HEf/A+k2[Ff/p] 



[P]t 



l^ATP 



[pATP 1 
l^ADP.Vil 



[P], 



k,C„ 



k,C„ 




ADP-V, 



ATP K,m ^s^^ATP 



Figure 2. Alternating cycles for Pgp. Scheme based on the original 
proposal of Senior et al. [25] that includes the coupling of two Elemental 
Cycles of ATP hydrolysis, the trapping reactions with Vj, the priming 
reactions with ATP, the priming reactions with ADP, the ADP- 
dependent Vj trapping reactions, and the interconnecting reactions 
between the ATP and ADP trapping pathways. E and F represent two 
ligand-bound isoforms of Pgp (P, the bare enzyme) with the ability to 
hydrolyze ATP in NBD1 (superscript position) and NBD2 (subscript 
position), respectively. The nomenclature for the rate constants 
corresponds to that defined for the Elemental Cycle (Figure 1, rate 
constants in Table 2). The cycle (shaded area) supplemented with the 
blue reactions corresponds to the PE Alternating Cycle (rate constants in 
Table 3). The addition of the red reactions defines the Extended 
Alternating Cycle (rate constants in Table 4). 
doi:10.1371/journal.pone.0098804.g002 



with P, E and T defined in Figure 2, and the vector k upgraded to 
include (kQ,k_Q). Furthermore, additional reaction paths were 
added to account for the Extended Alternating Cycle (Figure 2, 
including red reactions), for which the following new variables 
were defined: the fraction of single-nucleotide trapped species, Ti, 
and the fraction of two-nucleotide trapped species, Tji, given by 



Results 

Setting and Evaluation of the Elemental Catalytic Cycle 

The steady-state solutions of the biochemical variables for the 
Elemental Cycle correspond to the following expressions 



ri=T(k,Co,[P],)= 



r„=T(k,Co,[p],): 



-[Fa 



[Ea 



[P\, 



[ + [Fadp.v} 



[P\t 



k,C„ 



k,C„ 



.\ATP\ 



K,, 



, \ADP] [Pj] [ADP] 
-ADP Pj icADP 



K-1 



+ 1 + 



\ATP\ 



(1) 



by upgrading the vector k to include <fcia,^-iaj^4a>^-4a>^5a>^-5a)>- 
The analytical solutions were obtained using the computational 
algebra package GROEBNER included in Maple 15 (MapleSoft 
Inc., Waterloo ON, Canada), while the general-purpose simula- 
tion package SCoP 3.5 (Simulation Resources Inc., MI, USA) was 
used for numerical integration. For some plotting and fitting 
procedures, OriginPro 8 was used (OriginLab Corp. MA, USA). 

The solutions to these functions were used to kinetically evaluate 
the self-consistency of the models in accounting for aU the 
experimental evidence, using as a biochemical system hamster Pgp 
in the presence of the substrate verapamil, in either native plasma 
membrane, or purified and solubHized in detergent. Table 1 
presents a summary of the values reported for hamster Pgp from 
the main papers in the field, including (i) phenomenological 
parameters such as /C^o, K'^^ , K"'''' , T„, v„ax in steady-state 
kinetic experiments, or k„i,^ in time-course experiments, from both 
catalytic activity and trapping experiments; and (ii) thermody- 
namic and kinetic parameters such as K^i, K„,, K, and k^at- 



I [Vj] \ADP] 



\ATP\ 



{ADF\ \Pj\ {ADP\ \Vj\ {ADP\ 

KjDP + ZPJ KjDP ^ ^ K- 
d d \i 



ADP 



+ 1 + 



[Pi] 



lATP] 



(2) 



with the steady-state concentration vector defined by 
Css = <[ATP],[ADP],[Pi],[Vi]}. Eqs. 1-2 can be evaluated under 
the assumption that Css is almost identical to Co under the 
experimental conditions, e.g. using low monitoring the initial 
rate of activity, and/ or measuring the initial rate in the presence of 
an ATP-regenerating system. In both expressions, the thermody- 
namic parameters are defmed by 



(3) 
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and the kinetic parameters by 



(a) k„, = (b) K,„ = ^ ( Kf,^'' + ^"j (c) Tff '■ = ,9c. (d) <' = ,9o<.<' (4) 



where the factors a and (p are defined as 



(a).^|l(b).9../l + ^ + ^) (5) 



From these basic parameters we derived the phenomenological 
parameters, K"'''' and /C50, for the observed steady-state turnover 
rate and trapped fraction. However, in the absence of rapid kinetic 
data for Pgp-mediated ATP hydrolysis, it was necessary to use an 
arbitrary setting for the vector k. Although, the exact values of the 
vector k are unknown, the settings specified liere were designed to 
reproduce as closely as possible the pubhshed experimental data in 
Table 1 . The rationalization of the Elemental Cycle rate constants is 
as follows: 

(i) The setting of (p was a key element in the modelling. By 
definition, its value is >1, depending on the values of fc_3 
and fc_4 relative to (Eq. 5b). Since P; has low aflfmity 
(see (v)), fc_3 will be high (fc_3>>l s ), so the second 
term in Eq. 5b will be negligible. In addition, fc_4was set to 
the same value as kj (see Discussion). From both 
considerations, (p takes the value of ~2, however, its 
validation comes from the reciprocal constraints imposed 
by the interplay of the thermodynamic and kinetic 
parameters (Eq. 4). 

(ii) The observed catalytic rate, A:^^, SlOs^' is equivalent to 
fmaar of ~ 1 0 (imol Pi mg~'min~' [8]. Since two ATP are 
hydrolyzed per Pgp in the fuU catalytic cycle (see Alternating 
Cycle below) the rate constant for the hydrolytic step was set 
to k2 = 2G s"' (Eq. 4a). 

(iii) In the absence of products and inhibitors, Eq. 1 presents 
the characteristic hyperbolic behavior observed for the 
ATP dependence of ATP hydrolysis by Pgp, according to 

with parameters within the range reported: a maximal 
£ turnover rate of kcat, and a high consensus Michaelis- 

E _: Menten constant of .^^,,,=588 |iM for verapamU-stimulated 

!_ I Pgp [5,8]. The low affinity of ATP was set at 

^ s /s:^^^= lOOO^M, based on a reference value of ~1 mM 

for the effect of ATP on inhibition of labeling of NBD 1 by 
8-azido-ATP [16]. From the selected Km and AT^"" values, 
^ o 5 ^iid f = ~2, the ratio ki/kx is predicted to be high (~200, 

"g i 15 ^"^^ ^^)' from which the rate constant for the association of 

Is ATP is estimated to be A, =0.1 LiM"'s"'. 

^ ™ gi ^.A^ 8 (iv) ADP was reported to compete with ATP for the 

.2 ^ S 5 ^ § nucleotide-binding site [5,9] . Effectively, the mathematical 

§ .a o 1 -^^ model predicts pure competitive inhibition behavior of 

^''■S:^ ^ S 1 ADP on ATP hydrolysis (Figure SI) according to 



< ^- 

.C 4-" 

4-' (TJ 



^ > 



1 T s B k ,\ATP\ 

v,,= """^^X^ when[P,] = [V.]=0 (7) 



d 
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Table 2. Rate constants for the Elemental Catalytic Cycle. 





Process 


Parameter 






ATP association 






0.1 nM^'s"' 


ATP dissociation 






100 s"' 






1000 (iM 




ATP hydrolysis 






20 


P] association 


*3 




5 mM"^"' 


Pj dissociation 


k-. 




500 s"' 




<• 


100 mM 




ADP association 


(r4 




0.04 nM^'s"' 


ADP dissociation 






20 5-' 






500 fiM 




Vj association 






0.015 nM^'s"' 


Vj dissociation 


k-s 




0.020 s"' 






1.3 |iM 





Rate constants defining the vector kfor the reaction scheme shown in Figure 1. The nomenclature of the subscripts is follows: (±1) for the ATP equilibrium, (+2) for the 
hydrolytic step, (±3) for the Pj equilibrium, and (±4) for the ADP equilibrium. A positive sign is used for association reactions, a negative sign for dissociation reactions. 
doi:l 0.1 371/journal.pone.0098804.t002 



with K"PP increasing with ADP concentration. Eq. 7 states 
that the ADP inhibition constant for ATP hydrolysis, 
Kf'^'', indeed corresponds to the ADP affinity, K^^^ . 
Thus, for a given K„,, the ADP affinity is constrained by 
the observed K'^pp . Herein, ifj'^^was set at 500 ^M, which 
is compatible with the experimental value. 
Given Kf' = 204 mM, which is close to the reported value 
[14], Figure 3 shows the simulated output of the 
hydrolytic activity when [ATP] and [PJ were varied, 
according to 



Inhibition of ATPase activity following reaction with V; has 
been extensively studied. It was demonstrated early on that 
the trapped species is the long-lived Pgp=ADP=Vi 
complex, independent of the nucleotide used, and that 
the release of ADP correlates well with the slow 
reactivation of the enzyme [23]. The initial rate of ATPase 
activity, measured after rapid (~30 s) removal of unbound 
ligands, is approximately proportional to the relative 
concentration of untrapped enzyme. From Eq. 1, the 
ATP dependence of Vi inhibition is described by 



vss = 



1 + 



[ATP] 



when [ADP] = [Vi]=0 (8) 



vss = 



[ATP] 



-[ATP] 



when [ADP] = [Pi] =0 (9) 



' constrains the value of the affinity constant for binding 
of P; to the complex E^^^. Thus, given the value of ^ and 



According to Eq. 9, /Cj'Jj 



oc= 1 
4c). 



Ky is estimated to be 100 mM (Eq. following equation for [PJ 



corresponds to the solution of the 



Table 3. Rate constants for the priming reaction of the PE Alternating Cycle. 





Process 


Parameter 




k 


ATP association 


ko 




10 fiM^'s"' 


ATP dissociation 


k-o 




50 S-' 






5 \iM 





Rate constants defining the vector k in conjunction with the rates in Table 2, for the blue reactions In Figure 2. The nomenclature of the subscripts is as follows: (±0) 

for the ATP priming equilibrium. 

doi:10.1371/journal.pone.0098804.t003 
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Table 4. Complementary rate constants for the Extended Alternating Cycle. 





Process 


Parameter 






ADP association 






0.04 nM^'s"' 


ADP dissociation 












50 |iM 




Vj association 






2x10"= HM^'s"' 


Vj dissociation 






1x10"" s"' 




< 


5 ^lM 




ATP association 






1 xlO"^ |iM"'s"' 


ATP dissociation 






3S-' 






30 mM 





Rate constants defining the vector k in conjunction with the rates in Table 2 and Table 3, for the red reactions in Figure 2. The nomenclature of the subscripts is the 
same as in Table 2, with the addition of the suffix a to identify this pathway. 
doi:1 0.1 371 /journal.pone.0098804.t004 



[23] based on the kinetics observed for ADP trapping 
kMDP] .„„ ki\Vi] .no 7/ 

[ATP] following the route E ^^-^ E'^'"^ -U- ^adp v. 

latter kn value would yield a forward rate constant 
_ 1 kcat[ATP\ /c_4 = 6xl0"'s-' (given ifj'^^ = 500/iM) which is far too 



-[ATP] 2 K„j + [ATP] low to be compatible with the k^^^ of 10 s ' observed for 

ATP hydrolysis. Therefore, for our simulation, the reported 
^4 was ignored and we used instead the value imposed by the 
mutual interplay among all the parameters, 
resulting in (viii) It has been reported that Vi-induced trapping is completed in 

about 10 s (^o6.5~0.3-0.4 s~ ') by the pathway 



jy([ATP]\Vi\=Q) ' \[ATP\ J (10b) fi-^^^'^' (Figure 1), with 200 nM ATP and Vi [23]. Thus, for 

the given ki (and and k2, reaction with Vi becomes the 



which is the same as for the Vi concentration dependence of 

ATP hydrolysis. Setting Kj' = l.33[iM, yields ^,'''=2.72nM 

(Eq. 4d) and the IC^^ corresponds to 4.32 |iM at 1 mM ATP (see 
Figure 4A). A value close to 4.0 |a,M was reported for half- 
maximal inhibition of Pgp ATPase activity by V; under the same 
conditions [23]. Similarly, the trapped fraction variable was a query 
in our analysis. The steady-state concentration of trapped enzyme 
defined in Eq. 2 follows a hyperbolic curve as the ATP 
concentration increases, according to 

{^^[ATP] 
Tss= ^ 'k^ when [ADP] = [P.] =0 (11) 

fi) 

At saturating Vi concentration (e.g. 200 |XM), the IC^^^ is 
indeed the "K"Pf" term in Eq. 1 1, which for the given Kp value 
yields IC^ff = 7.9 ^M, close to the reported value of 9.0 [23] 
(see Figure 4B for the untrapped fraction). 

(vii) Given the value of k_^, and Kj'^'' = 500 |J.M, was 
estimated to be 0.04 |j.M 's ' (Eq. 3c). This does not agree 
with kn = 1.2xlO^^|iM^'s^' as suggested by Urbatsch et al. 




l4 



[ATP] (|jM) 

Figure 3. Effect of Pf on Pgp ATPase activity. Semi-log plot from 
the evaluation of I'sslk c„ with Css = <lATF],0,[F,]^,fl} for = 0 (red), 
50 mM (green), 200 mM (yellow) and 1000 mM (blue). Inset: double- 
reciprocal plot with ATP concentrations ranging upwards from 100 |xM. 
Values of k are given in Table 2. 
doi:10.1371/journal.pone.0098804.g003 
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rate-limiting step, with a calculated lower limit of about 
0.0015 |iM 's However, taking into account the other 
pathway for breakdown of the intermediate E'^^^ in the 
absence of ^ £, with A:_4 = 20 s"'), /fs was here 

set to 0.015 nM"'s"'. 

Table 2 presents the assigned values of the rate constants (using 
the arguments above) for the 9 unidirectional reactions defining 
the vector k, and the derived dissociation constants, Kj, for the 4 
bidirectional steps (Eq. 3). Once defined, k was kept constant for 
the rest of the simulations and validations. 

The effect of P; on ATPase activity was a key element in the 
validation of the model. It was previously reported that Pi behaves 
as a mixed-type inhibitor of ATP hydrolysis [23], where 200 mM 
Pi reduces the apparent by 50%, while the apparent K„ is 
reduced by just 20% [14]. In contrast, according to Eq. 8, K'VP / 
k'cai (the slope of the lines in the Lineweaver-Burk plot, see 
Figure 3 inset) is independent of inhibitor concentration, since 
KSt I K'^'' =kcatl K„„ so that Pi behaves instead as an uncompet- 
itive inhibitor, which is incongruent with the reported data. 

Another inconsistency between the output of the model and 
experimental data comes from trapping with ADP. It has been 
reported that 1 5 |J,M ADP produces half-maximal inhibition in the 
presence of 200 |xM V; [23]. According to Eq. 2, the ADP 
dependence of the trapping is defined by 



V- I [^DP\ 



a 



when [ATP] = [Pi] =0 (12) 



1 + 



- [ADP] 



where at saturating Vi, the '"K^app" term corresponds to ICfif (as 
with ATP dependence, see Eq. 1 1). Thus at 200 mM Vi, given the 
K^'^'' and Kj' values, the calculated /C^o'°^ is 3.31 ^M, which is 
5-fold lower than the reported value [23]. 

Figure 4B presents the simulated nucleotide dependence of the 
untrapped (free) fraction, 1-Tss, for both cases given by the model 
in Figure 1. Fitting of the synthetic data gave a HiU number of 
n = 1 for both ATP and ADP, which is expected for binding of just 
one nucleotide according to Eqs. 11 and 12. However, the 
behavior reported experimentally was a steeper concentration 
dependence for both ATP and ADP [23]. 

The ability of Pi to protect Pgp from Vi trapping was also tested 
using the model. It was reported that in the presence of 200 |J,M 
Vi, protection by 200 mM Pi is negligible at 1 mM ATP, but 
becomes significant at lower ATP concentrations [23]. This 
differential Pi protection effect depending on ATP concentration 
could not be reproduced by the model in Figure 1. The 
evaluation of Tss as a function of [ATP] and [Pi] , in the presence 
of 100 |J,M Vi, is plotted in Figure 5, which shows lines of similar 
slope, and ATP dependence opposite to that observed experi- 
mentally, i.e. the slopes decrease at lower ATP concentration. 
From Eq. 2, Kf' decreases with increasing [ATP] according to 



K. 



-K^ 1 + 



Km 



[ATP] K 



when [ADP] =0 (13) 
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Figure 4. V; interaction with nucleotides in the trapping of Pgp. 

(A) Semi-log plot of the Vj concentration dependence of the trapped 
enzyme fraction with 1000 nlVI ATP (blue symbols) or 1000 ADP (red 
symbols), from the evaluation of Tsslkc^ with Css = <1000,0,0,[K/]> 
and Css = <0,1000,0,[K,]>, respectively." (B) Semi-log plot of the 
nucleotide dependence of the fraction of free enzyme, with 200 jilVI 
Vj and either ATP (blue symbols) or ADP (red symbols), evaluating 
T'sslk.Q, with Css = <[^ri'] ,0,0,200) and <0,[^£)/'],0,200>, respectively. 
Lines are the best fits to the Hill equation with n = 1. Values of k are 
given in Table 2. 

doi:10.1371/journal.pone.0098804.g004 

Another discrepancy between the behavior of the model and 
experimental data comes from the interaction of Vi and Pi with the 
E'-^^-Pcomplex. From the evaluation of Tss at 200 nM ATP 
(Figure 6A) or ADP (Figure 6B) as a function of [1^ and [/^] , the 
competitive interaction reported for these two anions is evident, 
according to 



PLOS ONE I www.plosone.org 



8 



June 2014 | Volume 9 | Issue 6 | e98804 



Kinetic IVIodels for P-Glycoprotein ATP Hydrolysis 



K l=K '\ 1 



Kn, , [Pi] 



[ATP] Kf 



when[ADP]=0 (14a) 



+ ^ when [ADP] =0 (14b) 



However, fitting of the synthetic data in Figure 6 to an 
expression using an effective inhibition constant for P;, K^lff, 
according to 



K. 



d.app ' 



--K 



[Pi] 



when [ADP]=0 



(15b) 



yields values of Ki^^^ =5l.8 and 45.4 niM, for trapping with 

200 |J.M of ATP and ADP, respectively, half the reported values of 

100 and 70 mM, respectively, after correction for ionic strength 

[14]. The experimental values might be matched by increasing 
p. 

Kj', but then the capacity of Pi to inhibit hydrolytic activity would 
be affected (see (»)). 

Considering the time domain. Figure 7A shows the time- 
course of the overall activity and formation of the trapped species, 
for a pulse of ATP and V;. Thus, evaluating T with 200 |J,M 
[^7T]„ and [Vi\„ (keeping both constant), the numerical simulation 
mimics the rapid formation of the trapped species (within 1 0 s) and 
the high steady-state fraction trapped that was reported in the 
Kterature [23]. However, the output of the model clearly disagrees 
with the reported transient kinetics of dissociation of the Vi- 
trapped state. Experimentally, upon removal of unbound ligands, 
the observed slow dissociation has ^o/,j= 1.4x10 s ' 
(Ti/j = 87 min) [23], which correlates with the recovery of ATPase 
activity by the pathway E'^'^'"^' ^ fi-^^^ ^ E 

(Figure 1). The step E^^^ E ocean, at the rate of 

/t'_4 

^_4=20s^', which is compatible with the turnover rate, and 
rules out this reaction as the rate-limiting step for recovery of 



activity. On the other hand, the step E 



E^'^P has a 



rate constant k- 



_s=k5Ky =0.02s-\ which is 140-fold higher 
than the observed dissociation rate. Thus, in order to match the 
observed kinetics of ATPase recovery either (i) the dissociation 
constant kJ' must be much lower than 0.01 (J.M (see Figure 7B), 
a value which is incompatible with the observed Kp for trapping 
with ADP and ATP (see above), or (ii) the association constant ks 
must be much lower than 0.015 s~', which is incompatible with 
the fast formation of the trapped species (Figure 7A). 

The slow recovery of ATPase activity from the trapped species 
might be explained by the existence of several hidden transitions in 

— £ADP Vi -pj-jjj possibility was 



the overall reaction E^^^ 



tested by adding a step with a low dissociation rate constant (< 
0.001 s ) to explain the slow backward reaction to form E^^^. 
Effectively, the pathway for the trapping reaction was substituted 



l-f 
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Figure 5. Protection of Pgp from Vi trapping by Pi. Plot of the P; 
concentration dependence of the trapped enzyme fraction with 
100 JIM Vj and different ATP concentrations, from the evaluation of 
7'sslk,c„ with Css = <[-4rP],,,0,[i=,],100> for lATP],= ^000 i^M (red), 200 
(green), 100 (yellow) and 20 \iM (blue). Values of k are given in Table 2. 
doi:10.1371/journal.pone.0098804.g005 



by 



;5 £ADP V, 



^6,* - 



£ADP V, ^ which describes con- 



secutive equilibria with K^^ and .Sr^'clissociation constants, 
respectively. In order to include a slow backward step and shift 
the equilibrium toward the species on the right, the new forward 
rate constant k(, was set to 1 xlO ^ % ' and the backward rate 
constant k_(, to 1 xlO * s ' (yielding ^^' =0.1), with a concor- 
dant increase of the Vi association equilibrium constant K^^ to 
~\QKP' . In this way, it would be possible to explain the slow 
recovery of ATPase activity, while the change in overall affinity of 
V;, Kj' , would not significantly affect the Kp for trapping with 
ADP and ATP. However, inclusion of this additional step could 
still not explain the slow inhibition observed with ADP, by the 

pathway E ^-^^^ ^^bp k, -^hk issue will be 

considered further below. 



Construction of the Alternating Catalytic Cycles 

In this section, we evaluate the Alternating Catalytic Cycle proposed 
by Senior et al. [25]. In our adaptation of the model (shaded cycle. 
Figure 2), the two equivalent forms of the enzyme, E and F, 
correspond to states of the enzyme with similar energetic and/ or 
conformational states that differ only in the hydrolytic properties 
of their individual NBDs. This notation is necessary to distinguish 
between the two-nucleotide species, E'^jp^F'^jp, according to 
their NBD hydrolytic activity, i.e. the ii-form is capable of 
hydrolyzing only the ATP molecule bound at NBDl (but not at 
NBD2), and vice versa for the i^-form, thus moving the enzyme 
symmetrically between both states. Two different models can 
account for the E/ F forms; in both it is necessary to include ATP 
binding at each NBD of the bare enzyme as a first step (priming 
reaction) to get the initial intermediates of the cycle: (i) starting 
from the same conformer of the enzyme, P, recruitment of the 
NBDs to the nucleotide-bound state occurs randomly, with the 
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Figure 6. Effect of P; on the V; dependence of trapping. Semi-log 
plot of the V, concentration dependence of the untrapped enzyme 
fraction incubated with (A) 1000 nM ATP or (B) 1000 nM ADP, from the 
evaluation of Tss\k.c„ with Css = <1000,0,[^,]<„[ !'',]> and Css = 



<0,1000,[/',],„[K,]>, respectively, for [P,]c = 0 (red), 200 |iM (green), 
1000 p.M (yellow). Values of It are given in Table 2. 
doi:1 0.1 371 /journal.pone.0098804.g006 



and 



probability of occupancy given by the intrinsic afTmities of each 
E {where binding takes place at NBDl) 



NBD, so that P 

ATP 

or F 5 



^ATP (where binding take place at NBD2), or (ii) both 
conformers of the empty protein {E and F) co-exist, each exhibiting 
its own constitutive binding properties {E allows binding at NBD 1 , 
while F allows binding at NBD2); they may or may not be 
kinetically connected by the equilibrium E <- ^ F. The kinetics 
of ATP hydrolysis and Vi trapping are identical in both models. 
For the sake of simplicity, we decided to work with the first model, 
with the conservation of mass given by 



with ij for: none, ATP, ADP, ADP-Pj and ADP- Vi, according to 
Figure 2. In this model, the intermediates E'^^^ and Fjijp exhibit 
the same properties regardless of their origin, whether from the 
priming reaction or a later hydrolytic event. In this particular 
implementation of the Alternating Cycle, transformation between the 
two forms was achieved by exchange of ADP from/ to the two- 
nucleotide intermediates to/from the one-nucleotide intermedi- 

ADP ^ , ^4TP ADP ^4TP . , 

, respectively. 



ates: 



pADP 
^ATP 



f J7^TP pATP 
^ATP and P^op ^ ^ E 



However, it would be equivalent to assign the transformation to 
either the hydrolytic step (which looks reasonable) or the P; 
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Figure 7. Time-course of ATPase activity and formation of 
trapped Pgp. (A) Transient behavior of ATPase activity (red) and the 
fraction of trapped enzyme (blue), evaluating r|]iQ_ with 
Co = <,[ATF]j,,0,0,[Vi]^} at the indicated concentration pulses 'of ATP 
and Vj. (B) Time-course of the fraction of trapped Pgp according to V, 
affinity. Transient behavior of the fraction of trapped enzyme on 
incubation with ATP and V,, evaluating with 
Co = <.[ATP]^„0,0,[Vi]^} for pulses of 200 |ilVI ATP and V, "of 50 s 
duration (not shown). Each curve corresponds to Kj' = 3 |iiVl (blue), 
0.1 |iM (black) and 0.01 nM (red). Values of k are given in Table 2; 
[Pit = 0.25 nM. 

doi:1 0.1 371/journal.pone.0098804.g007 
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dissociation step, since in eitiier case, tiie kinetic behavior of the 
system is the same. To maintain symmetry, this step was assigned 
to the dissociation/association of ADP, which is the last hydrolysis 
product to leave the NBD. In addition, it was necessary to include 
the trapping reactions with Vi for each half-cycle 

ipADP ^ " 



pADP V, , pATP 
^ATP ^ADP 



pATP N 
^ADPVJ 



It is remarkable that neither of these two obvious steps (the 
priming and trapping reactions) has been depicted explicitiy in any 
reaction scheme that considers both hcdf-cycles simultaneously. 
The former was added later for first time by Urbatsch et al. [32], 
who considered that both NBDs binds ATP independentiy 
(priming reaction) and then come together (dimerization) to form 
the species with two bound ATP (although their concept was 
different from the one proposed here, see Discussion). We describe 
this new kinetic model, with both priming and trapping reactions 
(grey cycle plus blue reactions in Figure 2), as the Partial-Extended 
(PE) Alternating (Catal^^tic) Cycle. Any differences between the 
properties of the PE Alternating Cycle and a tandem repeat of the 
Elemental Cycle, can arise only from these additional reactions steps. 
Therefore, we were interested in evaluating the influence of the 
priming reactions in the ATP dependence of several observables. 
The steady-state solutions of the biochemical variables for the PE 
Alternating Cycle correspond to the following expressions 



kM\ATPf 



|F,-| lADFl 



[r,] [ADP\ 



^[ATP]+(^l + ^ + ^^ 



[ATpy 



(16) 



[ATP]-' 



(17) 



with the steady-state concentration vector defined by 
Css = <lATP],[ADP],[Pi],[Vi]}. For the evaluation of Eqs. 16 
and 17, it was assumed that Css S Co, as explained earlier. For this 
model, the velocity and trapping equations are no longer 
hyperbolic. For example, the ATP dependence of the turnover 
rate now follows a quadratic equation given by 



kca,[ATPY 

KI^PK,„+K4ATP] + [ATPf 



(18) 



when[ADP] = [Pi] = [Vi]=0 



where kcat and K„ are the parameters corresponding to the 
previous model (the Elemental Cycle), and the ATP affinity of the 
bare enzyme is defined by K^'' = k^olko. It IS important to note 

that, whatever the relative value of , Eq. 18 can be 

adequately fitted to the function 



vss-- 



kcuAATPf 
K,„ + [ATP] 



^ when[ADP] = [Pi] = [Vi]=0 (19) 



other hand, if we consider a much lower Kj^^ {K^^ < <Kj^^), 
n approaches 1 and the deviation from hyperbolic is negligible, 
and only observed at very low ATP concentrations. Because the 
majority of reports describe Pgp ATPase activity as Michaelian, 
we set the value of K^^^ in the |J,M range. This value also matched 

the low K,i value for the poorly-hydrolysable analog ATPyS [29] 
and other experimental evidence [33] explained by the model (see 
Discussion). Thus, simulating the PE Alternating Cycle with the 
parameters in Tables 2 and 3, the fitting that describes the ATP 
dependence of activity is an effective single K„, of 596 )j.M for 
n = 1 , a value very close to that obtained for the Elemental Cycle 
(Figure 8A). 

The interaction with ADP is now no longer one of simple 
competition (Figure 8B), and is described by 



vss = 



kca,[ATP\^ 



Ki^'K„+K^ h+\0\[ATP] + [ATP] 



(20) 



when[Pi] = [Vi]=0 

unless we include a pathway for the reaction P <^^^> E''^^^ and 
its equi\'alcnt for the _F-form (shown in Figure 2 in red, but not 
considered at this stage). Nevertheless, at high ATP concentration 
(e.g. [ATP] >100nM »K^''), the behavior is apparentiy 

competitive, as the literature indicates (Figure 8B, inset), since 
the interaction occurs mainly inside the catalytic cycle (where both 
nucleotides compete for the vacant site in _£'^"'and F^j.p), and the 
concentration of the bare enzyme, P, is negligible at that ATP 
concentration (see Figure 9). 

As expected, the observed properties with respect to Pi 
remained constant, with a effective inhibition constant, Rf' of 
~200 mM, since the relationship between the phosphate binding 
step and the hydrolytic step is conserved b(-t\vccn the Elemental 
Cycle (or tandem repeats of it) and the PE Alternating Cycle. However, 
the double-reciprocal plot of the ATP dependence of activity (not 
shown) has an upward curvature given by 



1 

vss 



{ 



\[ATP\ 
when [ADP] = [Vi; 



1 + ^ 



[ATP] 



(21) 



revealing that indeed the slopes are Pi-independent (as for the 
Elemental Cycle), but are now affected by the ATP concentration. 

Vi also behaves similarly in the PE Alternating Cycle, inhibiting 
ATPase activity at low concentrations. As indicated previously, the 
Alternating Cycle by itself cannot explain the cooperativity found in 
the nucleotide dependence of Vi trapping. This cooperative 
behaviour arises because of the priming reaction in the PE 
Alternating Cycle. From Eq. 17, producing synthetic data for the 
untrapped fraction, l-Tgg, with parameter values ofK^^ = 5 

and Kj' = 1.33 )j.M, (Figure 10), and performing an unweighted 
fitting according to 



with 1 and being an effective Michaelis-Menten constant. It 
is interesting to note that if Kj^^ = Kj^^ , the deviation from 
hyperbolic is appreciable only at high ATP concentrations. On the 



ss'- 



[ATP]" 



: + [ATP]" 



when [ADP] = [Pi] =0 (22) 
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1/[ATP] (1/mM) 

Figure 8. Steady-state simulation of tlie PE Alternating Cycle. (A) 

ATPase activity. Semi-log plot of ATP turnover rate (symbols) from the 
evaluation of I'sslk.Q. with Css = <[.4rP],0,0,0>. The line is the best fit 
to a hyperbolic equation. (B) Inhibition by ADP. Double-reciprocal plots 
for ATP turnover rate from the evaluation of v'sslk.c with 
C^s = <XATP\,[ADP]^fifly for [ADP], = Q (red), 250 (green), 
500 |iM (yellow/) and 1000 |iM (blue), with ATP concentration up to 
100 |j,M. Inset: Double-reciprocal plots with ATP concentrations ranging 
upwards from 100 |ilVl. Values of W are given in Tables 2 and 3. 
doi:1 0.1 371 /journal.pone.0098804.g008 

we obtained a Hill number of n= 1.21 and KfJ^j^j or ICfff^ of 
~20 |a.M. Unfortunately, there is no experimental data published 
for hamster Pgp to compare with the HiU number obtained by 
simulation. For reference, IC^^'' for the closely-related mouse Pgp 
was reported to be 18 |iM, with n= IJ [32]. Since Eqs. 16 and 17 
describe with more accuracy the turnover rate and trapped 
fraction, they should be used to set the value of K^', rather than 
Eqs. 1 and 2, as was done in the previous section. 



For the priming binding reaction with ATP, 

k„[ATP] .j.p „ , 
r h orr^jp, the association rate constant was set to 

100-fold the value for the corresponding rate inside the cycle, 

J7ATP _|J — ^ £atp (^^d the _F-form equivalent), which is 
/:o = 10 ^M~'s~'. Thus, the priming reaction would not limit 
the establishment of steady-state catalysis. In addition, this 
relatively high value for the priming association rate constant 
permits further decreases to allow our model to explain the 
observed impairment in trapping behavior in some systems 
[29,34]. 

However, some important experimental data still remain 
unexplained according to the PE Alternating Cycle: (i) the slow 
kinetics of Vi inhibition with ADP, (ii) the slow kinetics of 
reactivation of ATPase activity, and (iii) the stoichiometry of 1 : 1 
Pgpinucleotide in the trapped species, where ADP is trapped with 
Vi. Indeed, according to the kinetic reactions in the PE Alternating 
Cycle (Figure 2, grey cycle plus blue reactions only), the trapped 
species should contain both ATP and ADP [E^jp^' and Fjl'p.y), 
since there is no direct pathway to release ATP before Vi. 
Furthermore, according to this scheme, the bound ATP would be 
hydrolyzed when the enzyme re-enters the cycle upon Vi release. 

As pointed out above, there is the need to add plausible steps 
that account for the observed kinetics of trapping and release of 
both nucleotides. If we now incorporate the red reactions. 
Figure 2 outiines a minimal reaction pathway, including (i) 




[ATP] (mM) 



Figure 9. ATP dependence of several variables according to the 

PE Alternating Cycle. Semi-log plots of the steady-state ATP 
dependence of the normalized concentrations of (i) one-nucleotide 
species (red): [E^'^''] + [F^„] with Css([ATP]fififi) (ii) two-nucleotide 
species (brown): [Eatp] + [F-jfp] with CssilATP],Ofifi) (iii) bare 
enzyme (green): [P] with Css{lATP]fl,0,0) and the relative hydrolytic 
activity, by evaluating vsslkc„f°'' ("^) inhibition (pink): 

Css(lATP],mO,0,0)/CssilATPio,OfiO) (v) P; inhibition (yellow): 
Css(iATP'],0,200fl)/Css{[ATP]flflM) and the normalized trapped 
fraction, by evaluating 7'sslk,c,B^°'^ trapped species (blue): 
Cc;>;([-4rP],0,0,200) (vii) Pi protection of Vptrapping (black): 
Css([ATP],0,1000,2QQ) Concentration values for Css are given in nM 
except for Pj, which are in mM. Values of k are given in Tables 2 and 3. 
[Pit = 0.5 nM. 

doi:1 0.1 371/journal.pone.0098804.g009 
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adaptation of the basic alternating cycle proposed by Senior et al. 
[25] (grey cycle and trapping reaction with Vi), (ii) the priming 
reaction with ATP (blue reactions), and (iii) the priming reaction 
with ADP and a pathway for the release of ATP from the two- 
nucleotide trapped species (red reactions), by an independent 
pathway different from the catalytic reactions. The model for the 
Extended Alternating (Catalytic) Cycle was simulated only by numerical 
methods, using the rate constants in Tables 2, 3 and 4. 

Figure 1 1 A presents the time-course of the concentration of 
total trapped species during exposure to Vi with ATP or ADP. 
Here, there is a noticeable difference in the rate of accumulation of 
the trapped species (the rising phase) depending on the nucleotide 
used, which can account for the different observed rates of 
catalytic inhibition for the two nucleotides [23]. In the presence of 

V; and ATP, by the pathway P £"4rp^^ 

^ATP E^jp^' {and the _F-form equivalent), the rapid 

trapping observed sets a relatively high value for ^5. However, 
the fraction of trapped enzyme is now distributed between the one- 
and two-nucleotide species according to the equilibrium 

^ATP^' E^^^'^'{and the F-form equivalent) defined by 

Kj^^ (see below). On the other hand, for V; trapping in the 

presence of ADP, by the path p'''^^^ e^^^^ eADP Y, 

(and the i^'-form equivalent), the slower observed rate sets a lower 
value for the Vi association rate, The assignment of this latter 
step as being rate-limiting comes from the need to keep Kj^'^^ in 
the high \iM range, so that k4a = k-nal K^^^ should not be low. 

It has been reported that regardless of the nucleotide used, the 
trapped fraction corresponds to the ATP-free species, i?'^^^ ^' and 
^ADP Vi [23]> therefore the dissociation constant /T^^'was set to a 
high value (high mM range). With both ATP and ADP trapping, 
after removal of the ligands (equivalent to turning off the 




10' lO'^ 
[ATP] (|jM) 

Figure 10. Steady-state simulation of tKie PE Alternating Cycle. 

ATP dependence of trapping. Semi-log plot of the ATP concentration 
dependence of the untrapped enzyme fraction (red symbols) on 
incubation with 200 |j,M Vj, from the evaluation of rsslkc„with 
Css = <[^7'i'],0,0,200>. Blue line is the best fit to the Hill equation, 
with n = 1.21. Values of It are given in Tables 2 and 3. 
doi:1 0.1 371/journal.pone.0098804.g01 0 



nucleotide and V; pulses in our simulation), the model exhibits 
characteristic reactivation of ATPase activity. This is inversely 
proportional to the trapped fraction of enzyme (the decay phase), 
and takes place at the same rate regardless of the nucleotide used, 
by the common path E'^'^'"'^' ^ f-^^^ ^ P (and the F- 

t^-5ci l<-A<t 

form equivalent). This requires that when trapping with ATP, 

the small fraction of two-nucleotide species, E'^jp ^' and F^^p.f^ , 

^ADP- v'l 



which can decay by two possible routes E^'^^ ^' 



^ATP 



k^iE'^jp (and the _F-form equivalent), must drain mainly by 

the left-hand pathway; thus k-i.^ was set to a higher value than 
k-^. With this setting, by making Vi dissociation rate-limiting, slow 
ligand release from the common species E^^^ ^' at a time constant 
of ~ 1 /A_5a correlates well with ATPase recovery for trapping with 
both ATP and ADP. 

In addition, it was interesting to investigate the effect of an 
additional pulse of ATP during the recovery phase. According to 
the model, in the absence of Vi the main trapped species is 
depleted by two possible routes: e^dp ^ j^adp r, 

k-5a 

kui[ATP] ADP'V 

E^j.p ' (and the equivalent for the F-form). Since recovery 

of catalytic activity has the same slow kinetics whether or not ATP 
is present [23], then ki.^<k-^.^ (Figure IIB). This low rate 
constant for ATP association, ki.^, agrees with the high Kj^^{ior 
the given k-i.^, in turn constrained as mentioned above). All this 
assumes that binding of ATP is the rate-limiting step in the 

f J i7/(Z)P F, '''"I'^^^l 17 ADP. V, J7ADP 

forward recovery pathway, ii ' ^ E^jp E^tp^ 

k-5 

which is in concordance with the relatively high setting of ^-5, 
constrained by the observed fast trapping and a |J,M value for Kj' , 
since = 



Discussion 

Understanding the catalytic cycle of Pgp is essential to elucidate 
its transport mechanism. In spite of the efforts of several research 
groups over many years in providing good quality experimental 
data, no detailed kinetic analysis has yet been carried out. 
Consequendy, some puzzling features of the system still remain 
vmexplained, including: cooperativity of ATP hydrolysis at low 
ATP concentrations; mixed inhibition of ATPase activity by Pi; the 
steep concentration dependence observed for Vi trapping with 
ADP/ ATP; the kinetics of Vi release from the trapped species; the 
kinetics of Vi trapping with ADP; the relative /C50 values for Vi 
trapping using ATP/ ADP; protection from Vi-trapping by Pi- and 
the detection of one-nucleotide trapped species. In this work, we 
present a quantitative evaluation of the currently accepted models 
for ATP hydrolysis and Vi trapping, and assess their ability to 
explain the accumulated biochemical data. Using analytical and 
numerical methods, we evaluated the steady-state and the 
temporal behavior of the two main observable variables, the rate 
of ATP hydrolysis and the concentration of trapped enzyme. 
Thus, the basic reaction scheme for hydrolysis proposed by 
Urbatsch et al. [23], and its implementation in the Alternating 
Catalytic Cycle [25], were tested for their ability to reproduce the 
kinetic behavior of these variables. 

The success and applicability of this mode of analysis depends 
critically on the set of kinetic parameters (rate constants) 
employed. Since such kinetic data does not currendy exist, we 
established a coherent collection of rate constants that simulta- 
neously matched both steady-state and temporal courses of 
all phenomenological and known thermodynamic properties 
describing catalysis and Vi trapping. This self-consistent set of 
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Figure 1 1 . Time-course simulation of the Extended Alternating 
Cycle. (A) Time-course of Vj trapping. Transient behavior of the 
trapped fraction evaluating TIj^q with C„ = (\ATP\^^fifl\Vj\ 'y and 
Co = <0,[^£)/']^,,0,[K,]^,> at the indicated concentration pulses of Vi and 
ATP (100 s; blue) or V; and ADP (1000 s; red), respectively. (B) Time- 
course of decay of the trapped species in the presence of ATP. Transient 
behavior evaluating r|k,c.. with Co = <[^rf]^,0,0,[K,]p> at the indicat- 
ed concentration pulse of V, and ATP (100 s), and a second pulse of ATP 



(100 s) during the recovery phase, by setting k 
(blue), and 10 ^ 

in Tables 2, 3 and 4. [P], = 0.5 nM. 
doi:1 0.1 371 /journal.pone.0098804.g011 



_ OO"^ (red), 10" 
^iM^^s"' (black). The remaining values of k are given 



parameters was obtained using the reciprocal constraints tliat 
impose: (i) the parameters tliat describe ATPase activity, i.e. /c^„„ 
A^,„ and Hill number n; (ii) reference values of for nucleotides 
and Pi; (iii) the kinetics and phenomenological KJ IC^^q of products 
(ADP and P;) and inhibitors (V ;) for hydrolysis and/ or trapping; 
and (iv) the temporal course of Vi trapping and post-trapping 
recovery of ATPase activity (which is invaluable). It should be 
noted that some of these parameters are species-dependent. For 
example, T^i^^ for trapping with Vi using ADP for mouse Pgp 
(ABCBlb/Mdr3) is an order of magnitude slower than that for 
hamster Pgp [32]. In this regard. Table 1 compiles most of the 
parameters and observables reported for hamster Pgp (ABCBla/ 
Mdrl). 



The Steady-State Properties of the Elemental Cycle 

As shown in Results, the output of this model is in agreement 
with the basic properties exhibited by an isolated half-cycle of ATP 
hydrolysis with respect to ATP dependence and competition by 
ADP. Our set of rate constants reported: (i) a high Michaelis 
constant (K,,, = 600/.(M) which, in combination with the relatively 
slow catalytic rate (A:^^„ S 10s"'), results in a low effective 
brmolecular rate constant k^^jK,„ = l.6 x 10''M"'s"'; (ii) inhibi- 
tion of ATPase activity by ADP at sub-mM levels 
[Kf^P s 500/(M); (iii) inhibition of ATPase activity by P; at high 
mM levels (iSCf' S 200mM); (iv) inhibition of ATPase activity by Vi 
at |J,M levels {Kp = 3\iM); (vi) nucleotide dependence of trapping 
at |J.M levels. AU of these values are the same order of magnitude 
as those reported in the literature for verapamil-activated Pgp 
(Table 1). 

However, this model could not account for either the mixed- 
type inhibition exhibited by Pi, or for the observed ATP 
dependence of its protective effect on Vi trapping [14,23]. Analysis 
of the steady-state expression in this model (Eq. 1) revealed that 
K^if and k"^f can be described compactly according to 



[ADP].f([Pi],k) + K^ 



when [Vi]=0 



- WW) 



when [Vi]=0 



(23) 



(24) 



where f and g are functions of [P,] and the vector k. Thus, in the 
absence of ADP, the ratio between both parameters at any Pi 
concentration would be constant. However, in the presence of 
ADP in the reaction medium, the numerator of Eq. 23 is not 
reduced to X^n, so the slope of the double-reciprocal plot is 
dependent on inhibitor concentration, a characteristic of mixed- 
type inhibition, as reported by Urbatsch et al. [23]. However, the 
explanation for the inhibition they observed is highly unlikely to be 
ADP accumulation following hydrolysis, since Pgp has a low 
catalytic rate, and the ATP concentration was kept constant 
during the experiment by a regenerating system. 

Analysis of trapping with ATP/ADP uncovered another 
discrepancy between the output of the modeled Elemental Cycle 
and experimental evidence. According to Eqs. 11 and 12, at 
saturating Vi concentration the /C50 values of both nucleotides are 
defined by 



ICi 



ly-ADP 



1 + 



lVi\ 



Considering that (i) the numerators follows the relationship 



K.>K: 



ADP 



and I 



K^'is always > K^' , since 



Vi . 



■-aD.K^' (Eq. 



4d) and C!(i9 = 1 + - — - H — > 1 (Eq. 5) for any value of the rate 

constants, the model cannot reproduce the experimental observa- 
tion that IC^^^ > IC'^Q ^ior any Vi concentration. To match the 
reported data K^^"" would need to be >2Ki„. Additionally, the 
steeper concentration dependence reported experimentally [23] 
obviously reflects the binding of two nucleotides in the fuU catalytic 
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cycle, and is in contrast with the Hill number of 1 obtained from 
the Elemental Cycle. 

The relationship between P; and V; also revealed an additional 
element that makes the Elemental Cycle unsatisfactory. By simulta- 
neously setting primary properties, such as Xf' and .^r,''' from ATP 
hydrolysis and /C50 for trapping with both ATP and ADP, for a 
given set of other fundamental properties (especially K,,, and 
Kj^^^), it was impossible to mimic the reported relationship 
between these two oxoanions. The Kp values for trapping with 
ADP and ATP reported by the simulation were half of the values 
obtained experimentally. These values could not be matched 
without changing the other reported properties, that is by either (i) 
increasing the competitive capacity of P; (decreasing K^'); this 
change increases the P; dependence of ATPase activity (i.e. by 
decreasing Kf' for activity), or (ii) decreasing the competitive 
capacity of V; (increasing K^'); this change affects the IC^^ and 
IC^o'^ for trapping. 

The Temporal Behavior of the Elemental Cycle 

Severed considerations indicate that the observed slow kinetics of 
ADP trapping cannot be used to estimate the rate of ADP binding, 
as suggested by Urbatsch et al. [23], since this would yield = 
1.2x10""' |a.M~'s~'. Such a low value for k4, for the given 
-^rf^^ ~ 500 |iM, would make ADP dissociation the rate-limiting 
step for ATP hydrolysis, even if Kj^^ is as high as 1.5 M. In this 
regard, decreasing k-4 below 0.04 s ' (keeping the other rates 
constant) has a profound effect on the catalytic cycle, decreasing 
the turnover rate and to unacceptable values. There is now 
ample consensus that catalysis is rate-limited in a concerted way, 
that is to say, there is no particular limiting step [18,35]. This can 
be rationalized if fc_4 is quite similar to k2, as long as the P; 
dissociation rate is large (fc_3 >> 1 s '), a requirement that is 
fulfilled due to the low affinity of P; for E^^ (and F^dp). Thus, the 
steady-state turnover rate would be limited only for the steady- 
state [E*^^ and [^Litp], which are in turn dependent on [477]. 

The Temporal Behavior of the Alternating Cycle 

In the case of ADP binding, it is not possible to incorporate 
additional unimolecular steps into the Alternating Cycle (as previously 
suggested [23]) without either affecting the overall Kj^^ while 
preserving the effective forward rate, or affecting the overall 
forward rate while preserving the overall Kj'^^- Nevertheless, 
Urbatsch et al. [32] considered fast binding of ADP followed by 
slow isomerisation but, again, inside the normal ATPase pathway. 
Our proposal on this issue, incorporated in the Extended Alternating 
Cycle, came from considering an alternative pathway for ADP 
binding (see the red reactions in Figure 2) outside the regular 
hydrolysis pathway. Thus, for V; trapping, by either the fast 
pathway using ATP or the slower pathway using ADP, the final 
intermediates are the same, E^'^^ ^ and F^j^p y,. This is the case 

since for the ATP pathway, the equilibrium E^^^ ^' <^^^> 
^ATP ^' i^^^ F-form equivalent) is almost completely shifted 
toward the left (i.e. Aia[477]/A;-ia <<1). Moreover, when the 
pulse of Vi and ATP is off (i.e. when ligand is removed), a rapid 
transition takes place toward the one-nucleotide trapped species 
(with k-u = i s"'). Thus, E^'^'"^' and E'*^'"^' would be the 
trapped species present in the gel filtration column eluate, as 
observed experimentally [23]. 

Both one- and two-nucleotide trapped species have been 
captured for hamster Pgp, depending on the nucleotide analog 
and inhibitor employed. In the presence of fluoroaluminate 



(AIF4), Sankaran et al. [36] identified two nucleotides as 
trapped/bound (ADP/ ATP or 8-azido-ADP/8-azido-ATP) when 
Pgp was incubated with ATP or 8-azido-ATP. In contrast, only 
ADP/8-azido-ADP was trapped in the presence of V;. Beryllium 
fluoride (BeFJ behaves similarly to Vi in combination with 
nucleotides/8-azido-nucleo tides [37]. Thus, the geometry of the 
transition state, dictated by the divalent cation, the inhibitor and 
the nucleotide analog used, determines the properties of the NBDs 
and the interaction between them, (i.e. Kj^^). 

The inclusion of an independent pathway for ADP binding thus 
offers important advantages in explaining the cycle of catalysis and 
trapping. However, a complete reaction pathway should, in 
principle, consider sequential binding of two ADP molecules, as in 

P ■ ^'^^> E^'^^ <^'> E^^^ (and the T'-form equivalent), in the 
same way that the PE Alternating Cycle proposes binding of two ATP 
molecules. Indeed, Pgp can bind two ADP (one in each NBD) in 

the absence of ATP. Qu et al. [38] reported the binding of two 
TNP-ADP molecules to Pgp by fluorescence titration, and 
Tombline et al. [34] found a Hill number of 1.7 for ADP binding 
to Pgp catalytic mutants. In addition, the ADP dependence of Vi 
trapping displayed cooperativity [23]. However, the species E^^p 
and F^l^p can only be formed by incubation of Pgp with ADP 
alone, and thus they wiU not exist during the physiological catalytic 
cycle. Thus, the existence of distinct trapped species depending on 
the nucleotide used might account for the differential sensitivity to 
coUisional quenching observed for Pgp trapped with ATP 
compared to ADP [39]. 

In our simulation, the unusually small values assigned to the rate 
constants for the trapping pathways (i.e. A;±5a and ki.^; Table 4) 
should be noted: association rate constants for ligand-enzyme 
interactions are normally in the range 10"^-10^ HM"'s"\ 
However, these values were set in order to reproduce, within the 
minimal reaction scheme, the kinetic behavior exhibited during 
trapping and release experiments. For example, consider Vi 
release in the presence of ATP by the path 

j,ADP V, j,ADP V, ^ j,ADP. ^^^^^^^^ ^j^^^ ^j^^ 

ATP dissociation constant was effectively set to a high value, 
Kj^'' = 300mM. The setting of to a low value (the rate-limiting 
step) was due to the setting of k-^ to a high value, which was in 
turn based on kinetic analysis of the Elemental Cycle. As mentioned 
above, it is feasible to include additional steps in Vi release, 

^ATP ^' F^TPy t° allow assigning more reasonable values 

to these rate constants. The /:5a step, which explains slow trapping 
by ADP, can also be split into several conformational steps. Even 
ADP association with the bare enzyme can be slow due to the 
absence of constraints imposed by k^^^ in the regular catalytic 
pathway. In this regard, we found up to five transitions in TNP- 
ADP binding to Pgp under pseudo-first order conditions, with the 
observed time constants spanning 5 orders of magnitude, ranging 
from ms to tens of seconds [40] . 

The Concept of Alternating Catalysis 

Alternating catalysis, which was originally proposed by Senior et 
al. [25], integrated two Elemental Cycles in tandem. It arises because 
of a mutual interaction between the two Pgp halves that allows 
only one NBD to be active for a particular protc-in conformation 
(i.e. NBDl for E, NBD2 for E). Consequently, because binding of a 
second ATP to NBD 2 is required to enable NBDl to carry out 
hydrolysis (and vice versa), catalysis alternates between two 
Elemental Cycles. This characteristic is the crucial distinction 
between this mechanism and the Sequential Mechanism proposed 
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by Sauna and Ambudkar [30], where alternation of the two 
Elemental Cycles has its origin in the nature of the ATP binding step, 
rather than the hydrolytic step. Thus, for the latter model, the 
presence of ATP bound at a particular NBD is proposed to 
prevent binding of a second ATP at the other NBD. Biochemical 
and structural evidence supports the existence of a ternary Pgp 
complex with two nucleotides bound; the currentiy accepted 
model of catalysis is that each NBD carries out the catalytic cycle 
in turn, enabled by the complementary NBD with ATP bound. 

In the Alternating Cycle, during steady-state activity of the enzyme, 
at least one molecule of ATP is always bound (see Figure 2, grey 
cycle); each ATP molecule to be hydrolysed must bind to a 
previously formed Pgp-ATP complex. However, for a newly 
synthesized Pgp molecule in the cell, or at the beginning of an in 
vitro ATPase assay, the protein does not have any bound 

nucleotide. Thus, the priming reaction P <^^^> E^^^ (and its 
equivalent for the F-form) must necessarily occur. This step has 
possibly been ignored in the past because it is "obvious", but it is 
necessary' to include it explicitly to provide a pathway for the 
protein to enter the catalytic cycle. As discussed below, this 
additional binding reaction generates changes in the ATP 
dependence of any measured variable, and suggests some new 
concepts about the catalytic mechanism. The simple Elemental Cycle 
simulation obviously cannot report interaction in the nucleotide 
dependence of any variable for the intact Pgp, since only one 
nucleotide is involved in the cycle. However, for the Alternating 
Cycle, the observation of n> 1 for in vitro trapping with ATP arises 
because of nucleotide priming reactions (in the case of in vitro 
trapping with ADP, the simulation still reported n=l, since only 
one ADP binding event was considered in this case). When 
simulating the PE Alternating Cycle, the value of the Hill number 
obtained for trapping is dependent on the ratio between the two 
ATP afEnities and the type of coupling between the NBDs. Thus, 
for sequential binding of two ATP molecules, the Hill number 
ranges from: (i) n=2, when the catalytic sites present no binding 
interaction {Kj^^ = Kj^^) but show intcrdcj)<;ndence at the 
hydrolysis step (alternating catalysis, mutual exclusion of hydro- 
lytic activity); (ii) 1<«<2, for a negative binding interaction 
{K'^^^ > KjJ^) with again, inter-dependence of hydrolysis (e.g. 
alternating catalysis); and (iii) n-\, for mutual exclusion in the 
binding {K^^^ = co; i.e. after binding of the first ATP, binding of a 
second ATP cannot occur) and independent hydrolysis; which is 
the case for either uncoupled/isolated half-molecules (the Elemental 
Cycle) or the Sequential Mechanism (Elemental Cycles in tandem). For 
Pgp undergoing a complete catalytic cycle at both NBDs, as 
already discussed, option (iii) is discarded. However, due to the 
absence of any quantitative reports of the value of n for wild-type 
hamster Pgp, it is not possible to rule out either of the first two 
possibilities based on trapping experiments. However, several 
pieces of evidence point towards the second option (a negative 
binding interaction): 

(i) The Sequential Mechanism [30] proposed aUosteric control of 
the ATP binding affinities of the two NBDs. In this model, 
the alternating feature of the hydrolysis arises from the 
impossibility of a two-nucleotide species due to dramatic 
reduction in the binding affinity for a second nucleotide 
when one is already bound. Evidence was presented that 
correlated decreased affinity for drug with decreased affinity 
for nucleotide, thus, accounting for the release of both at the 
end of the catalytic cycle. Viewing this proposal using our 
model, the species 

^ATP 

must have low ADP affinity in the 
empty NBD (i.e. NBD2) to account for release of ADP from 



Fjl^ (see Figure 2). The low value of Kj^'' could be 

explained by extending that property to ATP. The Sequential 
Mechanism is, therefore, an extreme case of a negative 
binding interaction (interaction factor ^co). In our case, 
modeling the PE Alternating Cycle with an interaction factor of 
200 (if^"" = 2W)K^^) yielded n = 1 .25 for V; trapping with 
ATP and preserved the other properties of the catalysis and 
trapping. 

(ii) The observed Michaehs-Menten (n=l) behavior of the ATP 
dependence of hydrolysis requires either complete indepen- 
dence [K^^^ = 00 ; no binding of a second ATP) or a higher 
priming reaction affmity (K^J^ < <Kj^''). Although high 
affinity binding of ATP has not been reported, this is not 
conclusive, since the various reports are imprecise or 
incomplete. Plots of the nucleotide dependence of ATP 
hydrolysis by Pgp have often started from a relatively high 
nucleotide concentration (e.g. 50 |xM), thus missing details 
of the low concentration part of the curve. It should be 
noted that the inclusion of a high affinity priming reaction 
generates curves for ATP dependence that deviate only very 
slightiy from the single-binding model, so that it would only 
be perceptible in either log or log-log plots. In addition, the 
low concentration part of the curve could only be taken into 
account using a weighted fitting to a HiU model; a non- 
weighted simple Michaelis-Menten fitting would miss the 
high affmity component. 

(iii) An interesting report by Buxbaum [33], which measured 
hydrolysis of ATP in the nM range, reported signfficant 
deviation from hyperbolic behavior. Upward curvature in 
the log-log plot was observed at low ATP concentrations, 
with a breakpoint at ~10 |iM, which can only be explained 
by interaction between the NBDs during catalysis. In 
addition, the author reported that activation of ATP 
hydrolysis by verapamil occurred only at high ATP 
concentration, which might be reconciled with our model 
by adding priming cycle for ATP hydrolysis (i.e. hydrolysis of 
the one-nucleotide species) uncoupled from drug transport. 

The essential steps in the alternating mechanism proposed by 
Urbatsch et al. [32] are depicted in the cartoon in Figure 12A. 
The ATP binding reaction is conceived as a random process, 
producing the two-nucleotide intermediate (C) without any 
distinction in their binding affinities. Subsequentiy, this interme- 
diate chooses a pathway toward either Dn or Dc, depending on 
which NBD last hydrolyzed ATP. This model requires the 
intermediate C to lia\'c some type of "memory'", i.e. C must 
possess some intrinsic difference based on the last hydrolytic event, 
for example, a slight difference in the forward rate (C— *D) 
between NBDl (N-end) and NBD2 (C-end). However, by 
definition, C must be identical regardless of the branch used for 
the priming binding step, so that the next step would have to be 
randomly selected. This places Senior's Alternating Mechanism in an 
awkward position: in the forward step from C, there is no 
guarantee of alternation of the two half-cycles. 

In contrast, our proposal for the Alternating Cycle (Figure 12B) 
considers serjuential ATP binding with decreased affmity for the 
second nucleotide, to produce distinct two-nucleotide intermedi- 
ates, Cfj and Cq. In this model, alternation is guaranteed since 
there is no common intermediate; there is no need to propose the 
existence of memory for any species. The model in Figure 12B is 
equivalent to that shown in Figure 2, where one branch (blue) 
corresponds to the £-form of Pgp, and the other (red) to the in- 
form, and the intermediate A corresponds to the Pform. Thus, the 
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release of ADP and the transition between kinetics forms in 
Figure 2 {E — F), are represented by tlie transitions 
Fjsr-e — »-Bc and Fc-« — >-Bn in Figure 12B. 

In summar)', in our implementation of the Alternating Cycle in 
Figure 12B, the needed asymmetry for the alternation of the two 
paths is structural in origin (it is contained within the overall cycle), 
and arises from the reciprocal negative allosteric interaction 
between the NBDs. On the other hand, Senior's Alternating Cycle 
shown in Figure 12A is functional in origin (it is facilitated by the 
functioning of the cycle), since it depends on the "memory" of a 
particular intermediate' for the prc'\'i()us liydrolytic event. This 
feature ultimately arises because the model does not consider the 
priming reactions, and considers only the cycling part of the 
scheme. 

The Occluded State 

The concept of occlusion proposed by Tombhne and Senior [4 1] 
can be easily supported in our current model, as depicted in 
Figure 12B by the transitions C D. For this, the conforma- 
tional transition between the non-occluded (C) and occluded (D) two- 

nucleotide species would be represented by E^^p — i- E^^p in the 
kinetic schemes, where E denotes the non-occluded state and E the 
occluded state, with equilibrium constant A;. This transition is not a 
binding event, since there is no direct exchange (association or 
dissociation) of ATP; so that the apparent ATP affinity of the occluded 
species {K^"^''), would, in fact, be the overall ATP dissociation constant 
for the second nucleotide, as represented by the serial 

equiHbria, E^^^ ^'i^Z^d"' Ei^ Ei^P, witii Kf =KfP 

{\+Kt)^^. Thus, as occlusion progresses in the forward direction 
(A^>>1), the apparent ATP binding atfinity is significantly increased 
relative to the true "microscopic" binding afiGnity {K^^'' < Kj^^). 
However, this additional transitional step is not necessary to account 
for the experimental data reported with Pgp mutants and ATP 
analogs, as explained below. 

According to our interpretation, the work of Tombline et al. 
with Pgp mutants [34] might correspond to a pseudo-equilibrium 
binding titration of the bare enzyme, due to impairment in the 
hydrolytic rate constant, which reduced ki by a factor of 1000. 
Figure S2A shows the steady-state distribution at various ATP 
concentrations of the intermediates E'^^^ and Fatp, which closely 

matches the equilibrium P — ^i- ° E'^^^ 'i< — ' E^^p (and 
the i^-form equivalent). By decreasing both rate constants of the 
priming reaction (/cq and k^Q, keeping constant) the 

experimental data of Tombline et al. [34] could be simulated. 
After removal of free ligands, the reactions that describe the 

system, P ^ E^'^P ^ E^^p (and die F- 

k-o k-i 
form equivalent) predict occlusion of the nucleotide (equivalent to 
trapping without Vi). This arises mainly from formation of the 
species E^^^^ (and F^j.p) due to slower conversion of 
£ATP ^ p ^j^j ^jjQ increased formation by dissociation of 

two-nucleotide species E^pp »- ii'^"^^(and the F-form equiv- 

k-\ 

alent), which built up to a higher steady-state concentration 
because of greatiy reduced hydrolysis. We decreased ko and k-o 
by 1000-fold for the Pgp mutants and assumed that after passage 
through a gel filtration column (~30 s) almost all of the two- 
nucleotide species become one-nucleotide species (since 
A:_i = 100i"', ~2000A:_o)- Figure S2B shows the fraction of 
Pgp with retained nucleotide (i.e. occluded species) at various ATP 



concentrations, and reports an overall affinity for ATP (Kj^^) 
quite close to K^^^ , in the |tM range. In addition, the variable 
effect of several drugs on the steady-state stoichiometry of 
occlusion at a fixed nucleotide concentration [27] might be 
accounted for by a differential effect on K^^^ (see Figure S3). 

We recently reported that the binding to Pgp of the poorly- 
hydrolysable analog ATPyS exhibits a biphasic isotherm, with 

j^ATPyS ^ g ^jyj ^ATPyS ^ ^gj^ ^j^j^ highest 

affinity binding component assumed to correspond to the occluded 
state. However, our interpretation based on the current model is 

ATPyS 

that Kj^ may correspond to the affinity of the bare enzyme 
[KjI^p , see Figure S3). This would account for the Kf^^'^oi 
6 )j.M observed for the inhibition of the ATP hydrolysis by ATPyS 
[29], corresponding to competition with ATP for the high affinity 
site [Kj^P) of bare Pgp. 

In consequence, our model is compatible with the proposal of 
occlusion; but differs from it conceptually in the following way. 
The occluded state {E'^tp) identified experimentally has a tightiy 
bound ATP (|xM affinity) committed to hydrolysis, while a second 
molecule is bound to the complementary NBD. In our interpre- 
tation, this intermediate (Dn/g in Figure 12A) corresponds to the 
species E^jp and F^^p in Figure 8; it has one NBD with high 
affinity (e.g. NBDl and NBD2 for die E and F isoforms, 
respectively), and is represented by intermediates Cj^/q in 
Figure 12B. That is to say, in our model it is not necessary to 
include a conformational transition to the occluded state, since the 
high affinity site exhibited by this state actually corresponds to the 
site that bound the first nucleotide molecule to the bare enzyme. 
Consequendy, occlusion would not necessarily reflect an int t( ase 
in affinity of the NBD with ATP already bound, following binding 
of a second ATP in the c:omplementary NBD. Rather, it might 
correspond to the conformational change that enables the high 
affinity NBD to hydrolyze the committed nucleotide, thus 
preparing the enzyme for the hydrolytic step. This could occur 
concurrendy or after the binding of a second ATP, represented in 

Figure 12B as the transition Cj,j/c ^ Dn/c- The occluded 

state is easily incorporated into the kinetic scheme in Figure 2 as 

the transitions Ej^p <''^~> E^pp E^pp-^', however, as 

indicated previously, from a kinetic point of view it is not necessary 
to include this feature in our model. 

Conclusions 

The detailed analysis provided in this work underscores the fact 
that the mechanism underlying the kinetics of Pgp-mediated ATP 
hydrolysis must be much more complex than that proposed in 
previous models. Our goal was to incorporate the wealth of 
experimental data accumulated for hamster Pgp into a consistent 
kinetic simulation of the catalytic cycle. Implementation of the 
Elemental Cycle in the Alternating Mechanism (as originally proposed by 
Senior's group [25]) adequately explains (i) the time-domain and 
steady-state experimental data for ATP hydrolysis with respect to 
ATP, ADP and V; concentrations; (ii) the steady-state experimen- 
tal data for ATP/ ADP dependence of V; trapping; and (iii) the 
kinetics of V; trapping with ATP. However, it fails to satisfactorily 
explain (a) the effect of P; on ATPase activity; (b) the relationship 
between /C50 for ATP/ ADP on Vrtrapping; (c) cooperativity of 
ATP hydrolysis at low ATP concentrations; (d) the observed 
protective effect of P; on V;-trapping with respect to the /Cr,() for 
ATP/ ADP; (e) the steep concentration dependence observed for 
V; trapping with ADP/ ATP; (f) the kinetics observed for V; 
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Figure 12. Cartoons depicting the Alternating Cycle. (A) Random binding model adapted from Urbatsch et al. [32], (B) Sequential binding model 
proposed in this work (see Figure 2). White triangles represent ATP, black triangles represent ATP committed for hydrolysis; ADP-P, is shown in red, 
ADP in white. The subscripts of the intermediates (A to F) correspond to the N and C terminal halves of the protein. Closure of the NBD dimer 
interface is reflected in the fusion of both halves of the protein square. Flow through each half-cycle is represented by the blue and red arrows. 
doi:1 0.1 371 /journal.pone.0098804.g01 2 



trapping with ADP; (g) the kinetics observed for Vi release from 
the trapped-species; and (h) detection of species with only one 
trapped nucleotide. Development of the Extended Alternating Cycle 
allowed us to include additional kinetic steps to account for most of 
the deficiencies (c)-(h) of the original model (however, observations 
(a) and (b) still remain unexplained). Figure 9 summarizes the 
ATP dependence of several biochemical variables in the PE 
Alternating Cycle of Pgp, according to the parameters given in 
Tables 2 and 3. This proposed model introduces both priming 
and trapping reactions into the kinetic scheme, and is able to 
account for the observed high affinity of Pgp for ATP without any 
reference to the occluded state, thus avoiding assigning special 
properties to any intermediate in the cycle. A new interpretation of 
the occlusion phenomenon also emerges from the model. Future 
work will be needed to model a comprehensive reaction scheme to 
explain the complete data-set of biochemical observations. 

Supporting Information 

Figure SI Inhibitory effect of ADP on Pgp ATPase 
activity for the Elemental Cycle. (A) 3D plot from the 
evaluation of vsslk,Q, with Css = <[^rP],[^Z)P],0,0>, (B) Double- 
reciprocal plot from the evaluation vsslk,c„ with 



Css = <[^rP],[^Z)P]^.,0,0> for [ADP\, = Q (red), 250 (green), 500 
(yellow) and 1000 [iM (blue), with ATP concentrations ranging 
upwards from 10 ^iM. Values of k are given in Table 2. 
(TIF) 

Figure S2 ATP dependence of the concentration of 
several intermediates in the PE Alternating Cycle for 
catalytic mutants of Pgp. Based on Figure 2 and the values of 
k given in Tables 2 and 3, but substituting the following values: 
^0 = 0.01 nM"'s"', /l-_o = 0.05 s"' and ^2 = 0.02 s"' (a 1000-fold 
impairment in the original rate constants). (A) Concentration of 
intermediates: [P| (blue), + [F^jp] (red), and [Ei^p] + [Fj^] 

(green). (B) Fraction of Pgp with retained nucleotide, according to the 

function R^^\k,Q,M[E''^']HFATp]HEfT'p]HFiTp\)l[P], with 
Css = <[^rP],0,0,0>. The synthetic data from the model (blue 
symbols) were fitted to a HiU equation (red line), yielding 
r./, = 5.1 \xM and n= 1.01. 
(TIF) 

Figure S3 Effect of altering the affinity of the priming 
reaction on retention of nucleotides for catalytic mutants 
of Pgp. Kj^^ was altered by changing k-o, while imposing a 
1 000-fold impairment in the original rate constants (A^) = 0.0 1 |.lM~ ' s~ ' 



and kn = 0.02 s 



The fraction of intermediates with retained 
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nucleotide was evaluated by ^sslk,Cis defined in Figure S2) 
with Css = {[^ TP], 0,0,0), relative to the retained fraction for the 
original Kj^^^ = 5 |iM. Based on Figure 2 and the values of k 
given in Tables 2 and 3. 

(TIF) 

Figure S4 Simulation of the stoichiometry of trapped 
ATPyS based on the PE Alternating Cycle. ATP dependence 

of the stoichiometry of trapped nucleotide based on Figure 2, 
according to the function Qss\k,Qs = d^^^^^] + l^ATPys] + 
2[EiTpfs]+2[FiTPys])/lP], ^vithCss = <[ATPyS]M,Oy, for val- 
ues of k given in Tables 2 and 3, but considering 
^ = 0.01 nM"'s"', A:i = lxlO"*nM"'s"' and ^2 = 0.02 s"' (a 



1000-fold impairment in the original values). The synthetic data 
from the model (blue symbols) were fitted to a two-site binding 

[ATPyS] 

model (red line), yielding Qss~ 



[ATPyS\ 



5.Q2 + [ATPyS\ 



+ 



n99 + [ATPyS\ ' 
(TIF) 
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